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Introduction 

In this note we will consider the Potts ferromagnets and related systems on the square 
lattice. The Potts models are defined by the Hamiltonian 

^ = ~ ^ Jx,ySa x ,a y (1) 

with a x = 1,2, ... ,q and S ax!(Ty the usual Kronecker delta. Here J XtV is non-zero only when 
x and y are nearest neighbors and it is assumed that these couplings cannot be negative. 

Interest in the disordered version of these systems has recently been revived, in particular 
by J. Cardy and coworkers [Ca] who have discovered an apparent close connection between 
these problems and systems with random fields. For the disordered Potts models (among 
other systems) it was established in [AWi], [AW2] that for all temperatures, the energy is 
continuous. Thus, by conventional definitions, the magnetic ordering transition is continuous. 
However, as pointed out by Cardy in [Ca] - as well as in a number of public forums - for these 
systems it has not been established that all aspects of the transition meet the conventional 
criteria of a continuous transition: Vanishing of the order parameter, power law decay for 
correlations and infinite susceptibility. Here we establish that, at least at the self-dual points, 
these systems behave critically in the sense of all the above mentioned (with a power law bound 
for decay of correlations). Our results apply to a variety of systems under the hypothesis of 
a continuous energy density. 

The method is to employ graphical representations (e.g. the random cluster representation 
for the Potts models) and in fact applies to the non-integer cases provided the model is 
attractive (q > 1). In essence, the results here are complimentary to one recently proved 
in a paper coauthored by one of us [BC]. There it was shown that if (at some point) the 
energy density is not continuous then the discontinuity (a) is unique, (b) occurs precisely at 
the self-dual point and (c) coincides with the magnetic ordering transition. For these cases 
the picture is, by and large, complete. Unfortunately, for the continuous cases, our methods 
do not rule out the possibility that the self-dual point is simply a critical point in the interior 
of a critical phase. (Nor does it rule out the possibility that at the low-temperature edge 
of this purported phase, the magnetization exhibits a jump akin to the Thouless effect in 
one-dimensional long-range systems [T], [AY], [ACCN].) Nevertheless, taken together, the 
two sets of results imply that the self-dual point is always a point of non-analyticity. 

This work will be organized along the following lines: We will start with the uniform q > 1 
random cluster cases which are the simplest illustration of the basic method. Next we will 
treat certain straightforward generalizations, e.g. the Ashkin-Teller model and in the second 
section we treat the disordered Potts model. 
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Uniform Systems 

The random cluster model. We shall begin by setting notation. Consider the random 
cluster model on some finite connected AcZ 2 . If a; is a configuration of bonds, the probability 
of c<j, in the setup with "free" boundary conditions is given by 

« R N Mq c 'M (2) 

Where N(u>) is the number of "occupied" bonds of the configuration and C^(uj) is the number 
of connected components. If q is an integer, this is the representation of the of the model 
described by Equ.(l) - for free boundary conditions - with Jij = 1 and R = e 13 — 1. 

For other boundary conditions, the formula for the weights must be modified. Of primary 
interest are the wired boundary conditions which, back in the spin-system correspond to 
setting each spin on the boundary to the same value. Then the formula for the weights of 
configurations is the same as in Equ.(2) with C^(uj) replaced by C w {uj) where the latter counts 
all sites connected to the boundary as part of the same cluster. 

In this note we will restrict attention to the random cluster measures that are (weak, 
possibly subsequential) limits of random cluster measures defined in finite volume. Here the 
boundary conditions that we will consider - essentially those that are handed down from 
the spin-systems - are defined as follows: In volume A, the boundary, <9A is divided into k 
disjoint sets. All sites of the individual sets are identified as the same site. Thus the boundary 
consists of k effective sites (components) vi,...Vk- No interior connections between vi,...Vk 
are permitted. (I.e. any configuration with such a connection is assigned zero weight.) All 
interior sites connected to the same boundary component are considered as part of the same 
connected component. Finally, we will allow couplings between boundary sites and their 
neighboring interior sites to take arbitrary values in [0, oo). And, of course, we will also 
consider arbitrary superpositions of all of the above. 

For the case of free and wired boundary conditions, infinite volume limits, ergodicity etc. 
follow in a straightforward fashion from the monotonicity (FKG) properties of the q > 1 
random cluster measures. (In the disordered cases, some of these points must be rediscussed 
and we will do so at the appropriate time.) We will assume general familiarity on the part 
of the reader concerning these properties. Most of the relevant material can be found in 
[ACCN] or [BC]. However, if available at the time of reading, the authors highly recommend 
the forthcoming article [GHM]. 

The dual model, defined on the lattice A* that is dual to A, has weights of the same form as 
those in Equ.(2) with R replaced by R* = qR~ x . The general problem of boundary conditions 
for the dual model are a little intricate but for the purposes of this work, it is sufficient to 
note that the free and wired boundary conditions are exchanged under duality. 
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In these models (with integer q) the relationship between the bond density in the random 
cluster model and energy density in the spin-system is straightforward. In particular, let e x y 
denote the event of an occupied bond that connects the neighboring pair (x,y). Then, as 
shown in [CMi], 

{^(T x ,a v )%^ = ^ ( e z,J/) (3) 

where (— ) A '# denotes thermal expectation in the spin-system in boundary condition # = /, w 
(or, for that matter any other boundary condition). Thus, in these cases, continuity in the 
energy density is manifested as continuity in the bond density. Hereafter, we will focus on 
q > 1 random cluster models and use continuous bond density for our working hypotheses. 

Let us finally remark that for almost every R, the bond density is, in fact, a well defined 
concept. This point (which is fairly standard) has recently been detailed in [BC] so here we 
will be succinct. Consider the free energy, $(-R), defined here by 

$(R) = lim -|-logZ A# (4) 

A/"Z 2 |A| 

where Z A # is the sum of the weights in Equ.(2) and A /* Z 2 means a thermodynamic limit 
- a regular sequence of boxes. The function <&(R) is convex (as a function of logR) and 
hence has a left and right derivative for every R which agree for almost every R. At points 
of continuity of the derivative, there is uniqueness among the translation invariant states and 
the bond density in this state is given by the derivative of $ (with respect to log R) . At points 
of discontinuity, the upper value for the bond density is achieved in the wired state and the 
lower value in the free state. 

We are ready for our first result: 

Theorem 1. Consider the 2d random cluster models with parameter q > 1. Then if the 
self-dual point, R = R* = y/q is a point of continuity of the bond density, the percolation 
density vanishes. 

Proof. This result can in fact be obtained as a consequence of Theorem 2.1 in [BC]. For com- 
pleteness we will provide a direct proof. Here we will establish the contrapositive statement; 
i.e. assume that the percolation probability is positive at the self-dual point and show that 
this implies a discontinuity in the bond density. 

Percolation is defined in reference to the wired measures (and limits thereof). These 
measures are ergodic under Z 2 translations, respect the x, y-axis symmetry and have the 
FKG property. In short, these measures satisfy all the conditions of the theorem in [GKR] 
which forbids coexisting infinite clusters of the opposite type. Thus, with probability one, 
whenever there is percolation, all dual bonds reside in finite clusters. However, if there is 
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percolation (in the wired state) at the self-dual point, the same cannot be said in the limiting 
free boundary condition measure. Indeed, from the perspective of the dual bonds, this is a 
wired state. Hence, in this state, the dual bonds percolate and the regular bonds do not. It is 
thus evident that the limiting free and wired measures are distinct; ) is strictly below 

fJ^w(—)- By the (corollary to) Strassen's theorem (see [L] p. 75) this implies that the bond 
density in the wired state is strictly larger than that of the free state; the self-dual point is 
thus not a point of continuity for the derivative of $ □ 

Corollary. Under the hypotheses of Theorem 1, there is a unique limiting Gibbs state/ ran- 
dom cluster measure at the self-dual point. 

Proof. It was proved in [ACCN] (Theorem A. 2) that the absence of percolation (for q > 1) 
implies that a unique limiting random cluster measure and, for (integer q > 2) a unique Gibbs 
state in the corresponding spin-system. So this follows immediately. □ 

Let A and B denote disjoint sets in Z 2 . We will denote by {A < — > B} the event that some 
site in A is connected to some site in B by occupied bonds. Further, if D contains both A 
and B, we let {A < — > B} d denote the event that such a connection occurs by a path that lies 
entirely in D. For x and y (distinct) points in Z 2 , let g X:V = gf y (R, q) = ^ R ({x < — ► y}) be 
the probability that x and y belong to the same cluster. We will call this object the connectivity 
function. For integer q = 2, 3, . . . the connectivity function is equal to (or proportional to) the 
spin-spin correlation function. In these cases, the susceptibility and the average cluster size 
are also identified. For non-integer q, the geometric quantities are defined to be the objects of 
interest. These quantities are the subject of our next theorem which is a direct consequence 
of Theorem 1 . The proof below borrows heavily from the argument in [A] . 

Theorem 2. For self-dual q > 1 random cluster models with vanishing percolation probabil- 
ity, the function g XjV has a power law lower bound. Explicitly, if is the origin and L is the 
point (L,0) on the x axis then go,L, > \L~ 2 . Finally, the average cluster size is infinite. 

Proof. Since there is a unique limiting state then, in particular, the limiting free and wired 
measures coincide. Consider the L x L square centered at the origin which we denote by Sl- 
In every configuration, there is either a left-right crossing by regular bonds or a top-bottom 
crossing by dual bonds. By duality, in the limiting measure, these probabilities are both one 
half. 1 Letting Cl denote the left edge and TZl the right edge of the square, this implies that 

Yl 9x,y > \- (5) 

lr To ensure that this is strictly true, one must carefully construct the square so that it is exactly self-dual. 
However, for the arguments here, it is actually sufficient to observe that one of the probabilities must be 
greater or equal to one half. 
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Hence, for some (deterministic) x* E Cl and y* E TZl we have 

g x *, y * > 2^2- (6) 

This is, in essence the bound on the correlation function. For aesthetic purposes we will show 
that a similar bound holds for go, l but let us first attend to the susceptibility. Following the 
logic of Eqs.(5) and (6), there is an x** in Cl that is connected to ?Il by a path inside Sl 
with probability of order L _1 : 

^=v^ (K * m7 ^ }sJ > J_ (7) 

Regarding the point x** as being at the center of a square of side 2L and using translation 
invariance, we find 

^=^({0 — dS 2L }) > -L (8) 

i.e. X L = /j, q ' R= ^({0 < — > OSl}) > L~ l . This immediately implies a divergent susceptibil- 
ity/cluster size. Indeed writing 

* = Z>,« = E E so,* > E x ^ ( 9 ) 

x L xEdS L L 

this result follows. Finally, let us obtain our bound for the correlation function along the 
coordinate axes. Consider the event ({x** < — > TZl}s l along with along with its mirror image 
reflected along the midline of the square. I.e. a connection between Cl and y** by a path 
inside Sl where y** = x** + L. If these two connections occur in tandem with a top bottom 
crossing of 5x, we achieve a connection between y** and x** . These events are all positively 
correlated hence 

111 

#o,l = g y », X " > (2^21^2!^ ^ 10 ^ 

□ 

It is clear that the above generalizes to other random cluster systems. However since at 
present there are not too many examples of physically relevant models that satisfy all of the 
required conditions, we will be content with a small selection. 

The [r,s] -cubic (generalized Ashkin— Teller) model. Consider two copies of Z 2 with two 
sets of Potts spins: Tj E {1, . . . r} and Ki E {1, . . . s}. It is convenient to envision the model 
as two layers of Z 2 , the r-layer and the K-layer with the r-layer just above the /c-layer. In 
any case, the Hamiltonian is given by 

H = ~ E [^Ti.T-Ai,*,- + bS Ki>K . + c5 Ti7Tj \ (11) 
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where, as it turns out, we will be interested in the cases a, 6, c > 0. 

The dual relations for this model (at least for integer r and s) were derived some time 
ago in [dN], [DR] by algebraic methods. (Of course the special case r = s = 2 was derived 
much earlier starting, in fact, with [AT].) More recently, graphical representations for this 
model have been discovered [CMi] [PfV], (and see also [SS]) in which the duality is manifest. 
Consider bond configurations u> = (u> T ,u> K ) i.e. separate bond configurations in the r- and 
K-layers. As usual, we will start in finite volume. Let N(u T ) and N(u K ) denote the number 
of occupied bonds in the r- and k- layers respectively. Let N(u T V u K ) denote the number of 
edges where at least one of the r- or k- layers have occupied bonds and finally let N(uj t /\uj k ) 
denote the number of edges where both the r- and k- layers have occupied bonds. The 
graphical representation is defined by the weights 

W(oj) = A N ^ Ty ^^B N ^ TA ^^C [N ^ K) ~ N ^ T)] r c{ ^s c ^^ (12) 

where C(u T ) and C(u K ) are the number of connected components as in the usual random 
cluster problems. (And typically must be augmented with some boundary conditions.) The 
relationship between A, B and C and a, b and c is as follows: 

A= [(e /3& -l)(e /3c -l)]i (13a) 

e /3(a+b+c) _ e /3b _ e f3c _|_ ^ 



c = 



(e^ - 1) 



(13b) 



(13c) 



.(e^-l). 

In order for the graphical representation to make sense, we require 6, c > 0. However, this is 
not the case with a but it turns out that the FKG property - which we will need - only holds 
if B > A [BC] thus we actually require all couplings in Equ.(13) to be ferromagnetic. Under 
these conditions for the case r = s, b = c it was shown in [CMi] that there is a single ordering 
transition as the temperature is varied. 

The dual model is defined straightforwardly: edges of the dual lattice in, say, the r-layer 
that are traversal to occupied bonds are considered vacant dual bonds, those edges traversal to 
vacant bonds are the occupied dual bonds and similarly in the K-layer. The duality conditions 
are easily obtained from the weights in Equ.(12) (for the simple reason that V <-> A under 
duality) and the result is: 



A* = V" ( 14a ) 



(14b) 
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(14c) 



The analog of Theorems 1 and 2 for this system are readily established: 

Theorem 3. Consider the (r, s)-cubic model as described with A > B. Let $(A, B, C) 

denote the free energy similar to that defined in Equ.(4). Suppose that a self dual point: 
(A, B, C) = (A*, B*, C*), is a point of continuity for any first derivative of Then the per- 
colation probability in either layer vanishes. Let g T xy and X T denote the connectivity function 
and average size of clusters in the r-layer and similarly for g xy and X K . Then g^ L > \L~ 2 , 
X T is infinite and similarly for g£ and X K . 

Proof. It is convenient, but not essential, to restrict attention to the "plane" C = C*. Indeed, 
following the argument below, it can be shown that continuity with respect to A and B 
actually implies continuity with respect to C thus, for all intents and purposes, the (In- 
variable is out of the play. Continuity of the derivative with respect to B implies that at 
the point (A,B,C), the density of "doubly occupied" bonds is independent of state for any 
translation invariant state. Add to this continuity of the derivative with respect to A and 
(since N{uj t V uj k ) + N{uj t A ou K ) = N(lv t ) + N(oj k )) we may conclude that the total bond 
density is the same in every translation invariant state. However, we claim that this implies 
the same result for the separate densities. Indeed if the r-density were discontinuous, to keep 
the total density continuous would require a compensating discontinuity in the ^-density. But 
these densities are positively correlated; i.e. the discontinuities must go in the same direction. 
In particular, we could find a state (at the point (A,B,C)) where both densities achieved 
their lower value and another where they both obtain the upper value. 

Constancy of the bond density implies no percolation in either layer at a self-dual point 
which in turn implies unicity of the state. The rest of the argument follows mutatis mutandis 
the proof of Theorem 2. □ 

A loop related model. Our final example appeared in the context of loop models in [CPS]. 
Let u denote a bond configuration on Z 2 and let u denote the complimentary configuration 
on the dual lattice: If a bond of uj is occupied than so is the traversal bond and similarly 
for vacancies. (In other words, the vacant bonds of the dual configuration are the occupied 
bonds of the complimentary configuration.) The weights, in finite volume are given by 



We remark that from a technical perspective, unrelated boundary conditions for u> and u 
may be implemented. However it is natural to assert that if one is fully wired so is the other 



V{u) = L N ^s c ^s c ^\ 



(15) 
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and similarly with free boundary conditions. A derivation identical to the one for the usual 
random cluster model shows that the dual model is the same model with the parameter 

L* = j (16) 

along with the usual exchange of boundary conditions. For double-free or double-wired (as 
well as other) boundary conditions, the FKG property follows easily: 

Proposition 4. For free or wired boundary conditions, the random cluster models defined 
defined by the weights in Equ.(15) with s > 1 have the FKG property. Further, for these 
boundary conditions, if Ri > R2 and si < S2 the measure with parameters (Ri,si) FKG 
dominates the one with parameters (R2, ^2)- 

Proof. For the usual random cluster model, the FKG property follows from the FKG lattice 
condition [FKG]: Let 00 1 and 002 denote two bond configurations, uo\ A 002 the configuration 
of bonds occupied in both uo\ and 002 and u\ V 002 the configuration of bonds occupied in 
either. The lattice condition reads: fJ,(uo\ Au^M^i VW2) > h(loi) [1(002) which follows because 
C(oo\ A UJ2) + C(uji V U2) > C(uji) + 0(002) [ACCN]. In the present case, we need only apply 
this argument twice, once to C(uo) and once to C(Co). The FKG dominance follows by writing 
the one set of weights as an increasing function times the other. 

Remark. It would appear that the model under discussion is very close to the q = s 2 -state 
random cluster model. This follows by noting that for the former, C(oo) and C(Co) are iden- 
tically distributed. Then, we may write S <X<»)+C(G>) = s 2C(u) S c(u)-C{u) and suppose that the 

"fluctuations" C(Co) — C(uo) are (thermodynamically) small. However, at present there is no 
hard evidence of such an equivalence. On the other hand, the self-dual point can in fact be 
realized as the endpoint of the self-dual line of the symmetric (r = s, C = 1) cubic model 
corresponding to A — > 00 and B — > 0. Here the r-bonds may be taken to be the occupied 
bonds and the k's to be the vacants. The condition B = ensures that they cannot coincide 
while A = 00 implies one bond or the other actually is occupied. 

Theorem 5. For the model defined by the weights in Equ.(15), the results of Theorem 3 
apply; If R = s is a point of continuity for the bond density then the connectivity function has 
a power law lower bound and the average cluster size is infinite. 

Proof. Follows from the same arguments as the proofs of Theorems 1 and 2 □ 

Remark. In this model, the results of [BC] also apply: If there is any point of discontinuity 
for the bond density, that point must be the self dual point. (The cases of the Potts model 
and the cubic model for r = s were the explicit subject of [BC]; the identical arguments apply 
to the current case.) Thus, one way or another in all these systems the self-dual points are 
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points of "phase transitions". For large q, r, and/or s - at least in the integer cases - it 
is straightforward to show that discontinuities do occur. (Theorem IV. 2 in [CMj] covers all 
of these cases.) The difficulty is the opposite cases: establishing continuity. Indeed, to the 
authors' knowledge, the only non-trivial case where this has been done with complete rigor is 
the Ising magnet. However, the next section features systems where the required continuity 
has been guaranteed. 
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Quenched Potts Models 

For the remainder of this paper, we will deal exclusively with the g-state Potts model 
as defined by the Hamiltonian in Equ.(l); we will treat the case where the J XjV are (non- 
negative) independent random variables. (And ultimately to prove theorems along the lines 
of Theorems 1 and 2, we will need to focus on distributions that are self-dual.) We strongly 
suspect that with only minimal labor the forthcoming results could be extended to disordered 
versions of the various other models discussed in the previous section. But here we will focus 
on the minimal case. 

The approach in this work will be somewhat different from the usual mathematical studies 
of disordered systems: rather than looking at properties that are "typical" of configurations 
of couplings, we will construct, from the outset, the quenched measure - more precisely, 
the graphical representation thereof. When all the preliminaries are in place, this has the 
advantage of allowing a derivation that is essentially indistinguishable from the uniform cases. 
The disadvantage is that many of the "basic preliminaries" will require some attention. 

The quenched measure. Let AcZ 2 (or Z d for the duration of the preliminaries) denote a 
finite volume. In what follows, the inverse temperature (3 as well as the value of q (> 1) will be 
regarded as fixed and hence will be suppressed notationally. Let rj = {J x ,y \ (x, y) G A} denote 
a set of couplings. Let # denote a boundary condition on <9A. In general we will allow the 
boundary condition to depend on the realization of couplings so we will write ^(rf). (For the 
case of continuous variables J XjV we must also stipulate that #(?y) is a measurable function.) 
We let (— )A-#(r,) denote the finite volume Gibbs state (for this realization of couplings and this 
boundary condition.) Similarly, we may consider the random cluster measures l^\.^ rj - ) {—)- 
Our assumption about the J XjV is that they are i.i.d. non-negative variables. Let b(— ) denote 
the product measure for configurations of couplings and E{,(— ) the expectation with respect 
to this measure. Then the quenched measures are defined as the b-averages of the "thermal" 
averages according to (—) A .^^ and A t A-#(? 7 )( — )• Explicitly, if F(a Xl , . . - <T Xk ) is a function of 
spins (with x\,...Xk G A) then the quenched average of F is given by 

jF) A]# = E b ((F)l mri) ). (17a) 

Similarly, for a bond event A, 

JI A . i# (A)=E b ( f ,l #(r)) (A)). (17b) 

Most of our attention will be focused on the quenched random cluster measures as defined in 
Equ.(17b) - or the infinite volume limits thereof. Our first proposition establishes some FKG 
properties of these quenched measures: 
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Proposition 6. On finite A, let w and f denote the boundary conditions that are, respectively 
wired and free for all r\. Then the measures ~fi AjW (—) and ~fi A j(—) have the FKG property 
(in the sense of positive correlations). Furthermore, for any # = #(77), the wired measure 
~Pa,w(~) FKG-dominates the measure A*a,#( — ) 

Proof. For the FKG property, we will just do the wired case; the free case is identical. Let A 
and B denote increasing events (defined on the bonds of A). By the FKG properties of the 
measures At A;u ,(-), 

vlJAnB)>vlJA)vlJB). (18) 

Now we observe that for any increasing event C, the quantity n\. w (C) is an increasing function 
of 77. Indeed, if rji >- 772 (meaning J* > J% for each bond in A) then the random cluster 
measure with couplings 771 FKG dominates the one with couplings 772. But then 

Tl AjW (AnB) = ELIJAH B)) >M»l, w W A;w (B)) > 

J^))E b K ; J£)) = » A , w (A)tI A , w (B). (19) 

The same works for the free case (and various other boundary conditions that are independent 
of 77) . Finally, the stated FKG-domination is an obvious consequence of the fact that for each 
77, the measure ^ A)W (-) FKG-dominates the measure // A)#(7j) (-) □ 

As a corollary, we obtain the existence of infinite volume limits as in the usual random 
cluster cases: 

Corollary. The infinite volume limits Jif(-) and~p w (—) exist in the sense that if Ak / Z d is 
any sequence of boxes with Ajt+i D Ak then the limits ofJI Ak j(—) and~p AkW (—) exist and are 
independent of the sequence. These limiting measures are translation invariant and invariant 
under exchange of coordinate axes. 

Proof. The argument is exactly as in the standard proofs and is a consequence of the following 
observation: If Ai C A2 then for any fixed 77, the restriction of n A2 w (—) to Ai is FKG domi- 
nated by the wired measure in Ai. Thus the same statement holds for the quenched average 
of these measures. A similar sort of domination, but in the opposite direction is established 
for the free measures. The remainder of the proof is now identical to the derivations for 
uniform random cluster models (with occasional use of translation invariance and coordinate 
symmetry of b(— )). Such proofs have been written in many places (see e.g. [CMj] Theorem 
3.3) and need not be repeated here. □ 

We now demonstrate that absence of percolation is the correct criterion for uniqueness. 
Our working definition of percolation is fairly standard: 
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Definition. Let A C Z d be a finite set that contains the origin. We define 

-Poo = lim /J A (0 < — > dA). (20) 

A/"Z d 

We say there is percolation if P^ is not zero. 

Remark. It is obvious, by the considerations of the corollary to Proposition 6 that this limit 
exists. Further, if P^ vanishes, there is no percolation by any other criterion. Finally it is 
not difficult to show that P^ is exactly the spontaneous magnetization in the spin-system. 
Next we establish the quenched analog of Theorem A. 2 in [ACCN]. 

Proposition 9. // P^ = there is a unique limiting quenched random cluster measure and 
a unique limiting quenched Gibbs measure. 

Proof. Our proof will in essence be to show that any sequence of finite volume measures 
converges to the free measure. Let A denote any local increasing event. Let A denote a large 
(finite) box - the bonds of which determine the event A. Now consider a much larger box S 
along with some boundary condition #(ry); the measure ~p E ^(—) ma y be thought of as "well 
along the way" towards the construction of some infinite volume measure. 

Since the percolation probability is assumed to vanish, it is clear that if S is sufficiently 
large, then, for e > 

JI^(dA^dZ)<e 2 . (21) 

Thus if P A) s = {v I vl# (v) (dA <— > dZ) > e} then b(P A ,s) < c 

Now for any 77 G T>\^, with (il ^^-probability greater than 1 — e, there is a "ring" 
(separating surface) of vacant bonds in the region between dA and <9S. Conditioning to the 
"outermost" such ring gives us a measure which, in the interior of the ring, is equivalent to 
free boundary conditions on the ring. For any 77, this in turn is dominated by the measure 
with free boundary conditions on <9S and dominates (in A) the measure with free boundary 
conditions on dA. Thus, for 77 E T>a,e, 

(1 - e)[n\ f {A)] < < (1 - e)[»Z f (A)] + e (22) 

and hence 

(1 - e)[p AJ (A)] < JI^(A) < (1 - e)[p 3J (A)] + 2e (23) 

where the extra e comes from the e ^ Da,-- From Equ.(23) it is easy to see that all sequences 
of finite volume quenched measures converge to the limiting free measure. 

The argument for the uniqueness of the quenched Gibbs state follows from the above by 
noting that the thermal average of any local spin-function can be expressed as expectations 
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of random cluster functions (which themselves are finite combinations of increasing events.) 
This proves (a) the existence of a limiting (—) f (and for that matter a limiting (—) w ) and (b) 
that if the magnetization vanishes that this is the unique limiting state. □ 

The final result we will need is the ergodic property for the free and wired quenched random 
cluster measures. 

Theorem 10. The measures JI W (—) and]If(—) are ergodic under Z d translations. 

Proof. We will do the wired case, the free case is nearly identical. Let A and B denote local 
events assumed, without loss of generality, to be increasing. Let feZ d and let 7>(/3) denote 
the event B translated by f. We will show that lim^oo Ji w (A fl T?(B)) = Ji w {A)ji w {B) . 
By FKG and translation invariance we have, for any r, 

JlJAnT^B)) >Jt w (A)jI w (B). (24) 

Now consider \r\ large - far larger than the scale of the regions that determine the events 
A and B. Let s < \r\ be chosen so that A s , the box of side s centered at the origin and its 
translate by r, which we denote by TV(A S ), are disjoint but within a few lattice spacings of each 
other. Finally, let us consider an L that is very large compared with r; we will approximate 
JI w (AnT r (B)) by JI ALiW (AnT r (B)). By the FKG property, Ji AL W (Ar\T?(B)) is less than the 
corresponding probability given that all bonds on the outside of A s and T^(A S ) are occupied. 
But given these occupations, the measure inside A s is equivalent to wired boundary conditions 
on A s and similarly for T^(A S ). Now for each 77, the wirings make these interior measures 
independent. Thus we have 

»l L , w (AnT f (B)) < ^(^.(W (25) 

However as functions of 77, the two objects on the right of Equ.(25) are independent - they 
take place on disjoint sets. It is clear that A*t-(a ) W (^H^)) averages to Ji\ sjW (B) and thus 

»A L , w (AnT r -(B)) <JI AsjW (A)lI AstW (B). (26) 

Letting L — > 00 we get 

JI w (AnT?(B)) <JI As , w (A)jI As7W (B). (27) 

and hence 

Jim JI w (An T f (B)) < lim JI As , w (A)JI AsjW (B) = Jl w (A)Jl w (B). (28) 

This completes the proof for the wired case, the free case works the same way, here we use 
decreasing events for A and B. □ 
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Main results. We are ready for the disordered analogs of Theorems 1 and 2. However, in 
this case, we will not need to hypothesize the required continuity: this is the central subject 
of [AWi], [AW2]. Let us first briefly discuss duality in the disordered case. In the general 
setup, let 

Q 



J*(J) = log 



1 + 



(29) 



(I.e. e J — 1 = q/[e J * — 1]). Then what is needed, in the discrete case, to have j3 = 1 a point of 
self-duality is that b(J X:V = J) = b{J XyV = J* (J)). (A similar statement holds for continuous 
or other distributions) Indeed, if this is the case we see that the probability of bonds and 
dual bonds of equivalent strength are the same. Then, in finite volume, the probabilities of 
two coupling configurations that are equivalent under duality (including the usual exchange 
of boundary conditions) are equal. 

Sometimes it is convenient to parameterise the distribution and allow (3 to vary. For 
example, suppose there are two bond values J\ and J2 with b( J XiV = J±) = b(J XjV = J 2) = 1/2. 
Since temperature is back in the problem, we may assume, without loss of generality that 
Ji = 1 and write J2 = A with < A < 1. Then, in the A, (3 plane it is not hard to see that 
the model with parameters A, (3 is equivalent under duality to the one with parameters A*, 0* 
where 

Q 



(3* = log 



1 + 



and 



A* = 



log l + ^ti 



log 



1 + 



,A/3_! 



(30a) 



(30b) 



The system is self-dual (A = A*, = 0*) when (e p - l)(e A/3 - 1) = q. 

Theorem V . Consider a disordered Potts model of the type described at a self-dual point. 
Then there is a unique limiting state with zero magnetization. 

Proof. The proof is the same as the proof of Theorem 1 which we will recapitulate for conti- 
nuity. If the magnetization (percolation probability) were non-zero then in the free boundary 
state, the percolation probability for dual bonds would be non-vanishing. Proposition 6 and 
Theorem 10 allow us to use the result in [GKR]. Thus the free and wired states would be 
distinguished. But these are FKG ordered states so it would follows that the bond (and hence 
energy) density would differ in the two states implying a discontinuity in the energy density 
(or bond density). This however, is contradicted by the results of [AWi], [AW 2 ]. Proposition 
9 connects the absence of magnetization to uniqueness. □ 
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Theorem 2'. Under the hypotheses of Theorem 1' , in the limiting state the quenched corre- 
lation function satisfies 



/A 1 \> 11 



8L 2 



Further the quenched susceptibility, defined as 



q 



is infinite. 



Proof. Again this follows exactly the proof of Theorem 2 once we can conclude that the 
probability of a square crossing is one half. Although this follows from self-duality on "general 
principles" it is comforting to consider the square L x L square, Sl in the middle of a finite 
(but much larger) square with wired boundary conditions on the top and right and with 
free boundary conditions on the left and bottom. Then the quenched crossing probability is 
manifestly exactly one half and the thermodynamic limit can be taken, which gets us to our 
unique state, and we conclude that the probability in the limiting state is one half. □ 



L.C. takes great pleasure in thanking M. Aizenman; in general for many useful conversations over the years 
and in particular for some tips on this problem. And for the last reprint of [AW2]. 
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